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Abstract 



Building on a series of earlier papers fer-qc /9604007j ,9604008,9604009] , 
I investigate the various point-wise and averaged energy conditions in the 
Unruh vacuum. I consider the quantum stress-energy tensor correspond- 
ing to a conformally coupled massless scalar field, work in the test-field 
limit, restrict attention to the Schwarzschild geometry, and invoke a mix- 
ture of analytical and numerical techniques. I construct a semi-analytic 
model for the stress-energy tensor that globally reproduces all known nu- 
merical results to within 0.8%, and satisfies all known analytic features 
of the stress-energy tensor. I show that in the Unruh vacuum (1) all 
standard point-wise energy conditions are violated throughout the exte- 
rior region — all the way from spatial infinity down to the event horizon, 
and (2) the averaged null energy condition is violated on all outgoing ra- 
dial null geodesies. In a pair of appendices I indicate general strategy 
for constructing semi-analytic models for the stress-energy tensor in the 
Hartle-Hawking and Boulware states, and show that the Page approxi- 
mation is in a certain sense the minimal ansatz compatible with general 
properties of the stress-energy in the Hartle-Hawking state. 

PACS: 04.60,+v 04.70.Dy 



1 Introduction 

It is well-known that a quantum field theory constructed on a curved background 
spacetime experiences gravitationally induced vacuum polarization. This effect 
typically induces a non-zero vacuum expectation value for the stress-energy ten- 
sor [[j] ||, ||, ||, ||, @|. It is also well-known that this gravitationally induced vac- 
uum polarization will induce violations of at least some of the energy conditions 
of classical general relativity p[ f[ p[ p[ [To||. 
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In earlier publications, I addressed these issue in some depth for both the 
(3+1) Hartle-Hawking vacuum |8) and (3+1) Boulware vacuum ffl, and also 
presented an analysis for a (l+l)-dimcnsional toy model for which the exact 
analytic form of the stress-energy tensor is known p0| |. For the (3+1) Unruh 
vacuum there are additional subtleties that must be addressed. The complica- 
tions arise from two sources: 

1. For the (3+1) Unruh vacuum we do not have any approximate analytic 
expression for the stress-energy tensor. Nothing along the lines of the Page 
approximation for the (3+1) Hartle-Hawking vacuum or the Page- 
Brown-Ottewill approximation for the (3+1) Boulware vacuum |fL2[ , has 
yet been developed for the (3+1) Unruh vacuum. 

2. If we resort to numerical methods, using the numerical data of either 
Jensen-McLaughlin-Ottewill @ ||, |l| or Elster @, we must be very 
careful to adequately take into account a number of delicate and almost 
exact cancellations between various components of the stress-energy ten- 
sor. Even though individual components of the stress-energy tensor are 
often of asymptotic order 0[(2GM/r) 2 ], the particular combination rele- 
vant to the null energy condition [NEC] can easily be of asymptotic order 
0[(2GM/r) 3 ) or even smaller. 

To see where these cancellations come from, recall that the classic analysis 
presented by Christensen and Fulling Q is enough to tell us that asymptotically 
at large r 

(U\T^\U) = + 0[{2GM/rf]. (1) 

Note that I choose to work in a local-Lorentz basis attached to the fiducial 
static observers (FIDOS). The vector It = (1,1,0,0) is an outward-pointing 
radial null vector, L is the luminosity of the evaporating black hole, and the 
remaining terms in the stress-energy fall off at least as rapidly as 1/r 3 . 

If we look along the inward-pointing radial null vector [it = (1,-1,0,0), 
that is: £T = (-1,-1,0,0)], then 

ATj 

(U\T^\U) £7£7 = — + 0[{2GM/rf]. (2) 

This is enough to tell you that at large enough distances the Hawking flux will 
dominate, and that in this asymptotic region we will therefore satisfy the NEC 
along inward-pointing null geodesies. 

On the other hand, if we look along the outward-pointing null vector, then 
the term proportional to the luminosity (L) vanishes and so at the very least 

(U\T^\U) £pf = 0[(2GM/r) 3 }. (3) 
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I shall actually prove something considerably stronger. I shall prove that this 
quantity is negative everywhere outside the event horizon and that 



(U\T&\U) £pf = 0[{2GM/r) 5 }. (4) 

Thus for outward-pointing null geodesies the Hawking flux quietly cancels, and 
it is the sub-dominant pieces of the stress-energy that govern whether or not 
the NEC is satisfied. 

To get around these difficulties, I have acquired the Jensen-McLaughlin- 
Ottewill numerical data jl3|, |l4|, [ll| for the spin zero Unruh vacuum, and after 
suitable refinements (to be more fully described below), used this numerical 
data to construct a semi-analytic three-parameter model that globally fits all 
the known numeric data to better than 0.8% accuracy. Since the numerical 
data itself is not expected to have better than 1% accuracy, this is as good as 
we can reasonably expect, and further refinements to the model are not presently 
justifiable. 

Once I have developed this semi-analytic model, all discussions of the energy 
conditions in the Unruh vacuum will be phrased in terms of this model, and I 
will not further use the numeric data itself. 

As this work was nearing completion I became aware of a related though 
distinct analysis by Matyjasek jl7|. 1 compare and contrast my own analysis 
with that due to Matyjasek. 

Finally, in two appendices, I sketch a general strategy for constructing semi- 
analytic models for the Hartle-Hawking and Boulware states. 1 point out that 
there is a sense in which the Page approximation is the minimal ansatz com- 
patible with the known behaviour of the stress-energy in the Hartle-Hawking 
state, and show how to build a more general ansatz that does not disturb the 
low-order terms coming from the Page approximation. 

2 Vacuum polarization in Schwarzschild space- 
time 

2.1 Covariant conservation of stress-energy 

By spherical symmetry we know that for any s-wave quantum state \ip) 
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Here p, t, /, andp are functions of r, M, and %. Note that 1 have now set G = 1, 
and continue to work in a local-Lorentz basis attached to the fiducial static 
observers (FIDOS). I shall work outside the horizon and will make extensive 
use of the dimensionless variable z = 2M/r. 
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In the classic paper by Christensen and Fulling Q , it was shown that (by us- 
ing the equations of covariant conservation) the stress-energy tensor can, in the 
Schwarzschild spacetime, be decomposed into four separately conserved quan- 
tities. These four conserved tensors depend on two general functions of r and 
two integration constants. I choose to use a slightly different basis for this 
decomposition, and write the stress-energy tensor as 



<V|TH^> = [T trace r + [Tpressuref" + [T+f + [T-]*" '• 



Here 



[T t1 



w u = 



-T(z) + j^H{z) 

+ ^H(z) 



with 









0, 



(6) 
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(8) 



The conserved tensor [T tr ace]^' 1 ' depends only on the trace T of the total stress- 
energy tensor. Furthermore its trace is equal to that of the full stress-energy 
tensor itself: [T tr . ace ] M /i = T. 
Next, I define 



pressure] 



+2p(z) + ^G(z) 

+^ G ( Z ) 

+p(z) 

+p(z) 



(9) 



with 



G(z) = 



4 - 4 ) P{z) dz. 



(10) 



The traceless conserved tensor [T pressure ] tJ ' 1 ' depends only on the transverse 
pressure p(z) of the total stress-energy tensor. 

This decomposition makes sense if the integrals G(z) and H(z) converge, 
which requires mild integrability constraints on T(z) and p(z) at the horizon. 
These constraints are certainly satisfied for the Unruh (and Hartle-Hawking) 
vacuum state for which T\ z= i and p\ z =i are actually finite [|l8| 0. Thus 
we can write 



H(z) = l -T\ z=1 (l-z) + 0\{l-zf] 



G{z) 



-p\z=i {l-z) + 0[{l-zf 



(11) 
(12) 
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This tells us that near the horizon 



[Ttracer(z) 
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(13) 



and 



pressure 



+P\z=l 
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+P\z 



0(1 -z). (14) 



This is enough to imply that these two tensors arc individually regular at both 
the past and future horizons. (For the Boulware vacuum minor changes in the 
formalism are required. See Appendix B.) 

The function G(z) can be somewhat rearranged by an integration by parts. 
Suppose we define 



F{z) = 

Then it is easy to show that 



z 2 (1-z) 



dz 



p{z) 



dz. 



G{z) = - l —^p(z)-F{z). 



(15) 



(16) 



Near the horizon 



F{z) = 0[{l-zf]. 



(17) 



Doing this again requires only the mild assumption that p{z) is well-behaved at 
the horizon z — 1. Subject to this caveat the tensor [T pressure ] can be rewritten 
in the somewhat more convenient form 



Pi 



pressure] 



+p(z) £-F(z) 







-p{z) - £r x F{z) 






+p(z) 









+P(z) 



(18) 

In view of the definition of F(z) this particular form of the tensor makes it clear 
that a transverse pressure that falls of as 0(z 4 ) will have special properties. 
Finally, I take 
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[T±f = f± 44 = f± 



+ 1 ±1 

±1+10 







(19) 



These two traceless conserved tensors correspond to out-going and in-going null 
fluxes respectively. Furthermore /+ and /_ are two constants with the dimen- 
sions of energy density that determine the overall flux. In terms of the overall 
flux f(z) we have 



/(*) = (/+ - /-) 



1 



(20) 



Note that [T + ] is singular on the future horizon H + , and regular on the past 
horizon H~ . Conversely, [T_] is singular on the past horizon H~ , and regular 
on the future horizon H + . 

It is easy to see that this decomposition is equivalent to that of Christensen 
and Fulling. Perhaps the best starting point is equation (2.5) on page 2090 of 
reference [[18). The functions H{z) and G{z) defined above are not identical to 
those of Christensen and Fulling, but are linear combinations of the quantities 
appearing therein. Finally, the constants K and Q of that paper are in my basis 
given by K = (/+ — /_) 4M 4 , and Q = (—2/+) 4M 4 , respectively. I have chosen 
this particular basis because of its elegance, symmetry, and the ease with which 
I can adapt it to my semi-analytic model to be introduced later in this paper. 

Everything done so far works for an arbitrary quantum field in an arbitrary 
(spherically symmetric s-wave) quantum state in the Schwarzschild geometry. 
(In particular it holds for both the Unruh and Hartle-Hawking vacuum states, 
and with only minor modifications, also for the Boulware vacuum state.) I will 
now start to particularize the discussion. 



2.2 Conformally coupled quantum fields 

For any conformally coupled field, still in any arbitrary (spherically symmetric) 
quantum state, the trace of the stress tensor is known exactly and is given by 
the conformal anomaly. In a Schwarzschild spacetime 

T{z) = HPooZ & - (21) 

Here I have defined a constant (which in the Hartle-Hawking vacuum can be 
interpreted as the pressure at spatial infinity) by 

Po ° = 90(16tt) 2 (2M) 4 ' ^ 

The number £ depends on the particular quantum field under consideration. In 
particular, for a conformally coupled scalar field £ = +96. For other conformally 
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Table I. Anomaly coefficients. 
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All fermions are non-chiral. 



coupled quantum fields, just change the +96 to the appropriate magic number. 
The relevant coefficients can be easily deduced from the table on page 180 of 
Birrell-Davies || and are they are given in Table I. 

From the general definition (||) above, the function H(z) is easily evaluated 



H(z) = JL Poc (l _ z&) = I Poc (i _ z )(l + z + z 2 + z 3 + z 4 ). (23) 

Thus the trace contribution to the stress-energy is known exactly and analyti- 
cally, being given by a simple polynomial: 



[T t , 
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1 + z + z 1 + z 3 






9z 4 

1 + z + z 2 + z 3 + z 4 












(24) 

In particular this implies that once we have a semi-analytic model for the remain- 
ing free quantities p(z), /+, and /_, no matter how derived, we automatically 
have a model for the entire stress-energy. 



2.3 Unruh vacuum 

In the Unruh vacuum, we have a lot of additional information available beyond 
the rather general considerations given above |l8[ . First, since the Unruh vac- 
uum is to be regular on the future horizon we must have / + = 0, though it is 
permissible to have /_ ^ 0. Naively, this appears to forbid an outgoing flux, 
until we realize that there is nothing to stop /_ from being negative. In fact, 
this is exactly what happens, and so I shall define a positive quantity fo and set 
/_ = — fo- This is why we often hear the assertion that (at the event horizon) 
the Hawking radiation corresponds to an inward flux of negative energy. 
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Collecting the terms in the general decomposition (||) for the stress energy 
tensor, we have a this stage 



z 2 £ 

p(z) = -fo- + ±- Poo z 2 (l + z + z 2 + z 3 -9z 4 

1 — z 10 



+P(*) - — F{Z) - (25) 
z 2 £ 

r(z) = +/o— ^rPco z 2 (l + z + z 2 + z 3 + z 4 ) 

1 — z 10 

+pC0 + y^ f ( z )- ( 26 ) 

/(*) = +/o I ^- (27) 

At asymptotic spatial infinity we want the stress-energy to look like that of 
an outgoing flux of positive radiation [ fl8| . That is: we need to have p(z) — > f(z) 
asymptotically as z — > 0. From the Christensen-Fulling analysis, we know that 
asymptotically the transverse pressure goes as 1/r 4 fO(z 4 )]. So, picking off the 
dominant terms [0(z 2 )] at large distances we see that 

^ Poo - F(0) - /o = / . (28) 



That is 



Explicitly 



/o^-^. (29) 




I shall have occasion to use this result as a critical internal test for my semi- 
analytic model. 

Substituting this result for fo back into the general expression for the stress- 
energy, the various components are seen to be 
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^ Poc z 2 (1 + z + z 2 + z 3 



9z 4 
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-F(z). 



(32) 
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T[z) = + »oo Poo Z [1 + Z + Z + 

w V 20 2 ) 1-z 10 V 



z 3 + z 4 ) 



^ (Z) + I^ F(Z) - (33) 



2.4 Total luminosity: preliminary analysis 

The total luminosity of the black hole is now easily evaluated 

L = 4tt/o(2M) 2 . (35) 

More explicitly 



L = 4tt(2M) 2 



|j P°° + ^ / ( ^ ~ ~ ) d '- 



(36) 



This should be compared to the geometric optics approximation for the lumi- 
nosity, obtained from Stefan's law: 

L geometric— optics — ~Z & T S. (^7) 

The factor 1/2 comes from the fact that there is only one polarization state 
for a scalar field, while Stefan's constant is, in the current geometrodynamic 
units, simply a = tt 2 /60. [For SI purists: a = (7r 2 /60)(fc 4 /?i 3 c 2 ).] For the 
temperature we take the Hawking temperature T = Th = l/(8nM). [For SI 
purists: kTn — (he 2 ) /(8ttGM).] Finally, the effective radiating surface area 
is S = 4tt(3v / 3M) 2 = 108ttM 2 . [For SI purists: S = 108tt(GM/c 2 ) 2 .] The 
effective surface area is deduced from the well-known textbook result that any 
null geodesic of impact parameter less than a = 3^/3M will fall into the black 
hole (l^, ||^, . Equivalently, the geometric optics absorption cross section for 
null particles is cr C ross- section — 27-kM 2 Q. Pulling this all together 



L geometric-optics = SllT M 2 Poo = = 6.99411 X lO^M' 2 . (38) 

4UyoU7rM 

Restated as an estimate of the flux, we have 

81 

f geometric— optics — y^T Poo ■ (*^) 



2.5 Pointwise energy conditions: preliminary analysis 

From the analytic discussion presented so far, we have enough information to 
make a preliminary analysis of the pointwise energy conditions. I shall extract a 
very simple and robust result using a minimum amount of analytic information. 

We have already seen in the introduction that if we look along the outward- 
pointing radial null direction things simplify dramatically. So let's use what we 
have so far and calculate 



(U\T^\U) £fef = p-r-2/, (40) 

= + |pooZ 2 (l + z + z 2 + z i -Az A ) 

1 — z 5 

~2^-F{z), (41) 
z 2 £ 

= ~4foz + |p oc z 2 (l-z)(l + 2z + 3z 2 + 4z 3 ) 

1 — z 5 

-2-^F(z). (42) 
1 — z 

If we now work close to the horizon, and use the near-horizon asymptotic form 
for F(z) = 0[(1 - z) 2 } as deduced in eq. (|l|), then 



(U\T^\U)£+4 = -_*°- + 8/ + O(l-z). (43) 

So, provided only that the luminosity is positive, we will definitely have vio- 
lations of the null energy condition sufficiently near the event horizon. This 
automatically implies that the weak, strong, and dominant energy conditions 
will also be violated sufficiently close to the event horizon. 

Notice how minimal the input to this result has been: If the black hole 
radiates, and the conserved stress-energy has the correct asymptotic behaviour 
for the Unruh vacuum (both near the horizon and at spatial infinity), then the 
null energy condition must be violated sufficiently near the event horizon. 

This is in complete agreement with the celebrated area increase theorem 
of black hole dynamics p3|. From the area increase theorem we know that 
if a black hole radiates, and if that radiation back-reacts on the black hole so 
as to decrease its total mass, then the null energy condition must be violated 
somewhere in the spacetime. The area increase theorem is of course derived 
using totally different techniques. 

Naturally, once I have finished setting up my semi-analytic model, I will 
be saying a whole lot more about violations of the various energy conditions. 
In fact, I shall show that these violations of the point- wise energy conditions 
persist all the way from the horizon out to spatial infinity. 
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2.6 Model building: preliminary analysis 

Now we know, both from the analytical analysis of Christensen-Fulling |18fl , and 
the numerical analyses of Jensen-McLaughlin-Ottewill jl3| and Elster fllGfl , that 
the transverse pressure is well-behaved all the way down to the horizon. This 
suggests that the transverse pressure might be usefully modelled by some simple 
convergent power series in z = 2M/r. We also know that at large distances the 
transverse pressure must fall as z 4 . Suppose we normalize to p^, introduce a 
denumerably infinite set of dimensionless coefficients k n , and define 



p(z) = poo ^ k n z n . 

n=4 

Then, from the definition Jl5| ) of F(z) 



(44) 



F(z) = Poo^2(n - 4)k n 



In particular 



Thus 



1 - z n - 2 1 - z" 



n — 2 n — 1 



■E 

in — 3 



1 - z r 



[(to - 2)fc m+2 - (m - 3)A: m+ i] 



F(0)= Poc J2 



n— 5 



(n - A)k n 
(n-2)(n- 1)' 



/ £ ^ (n - 4)fc„ 

JO = Poo I ™ ~ 2^ 



20 



2(ra - 2)(n - 1) 



A brief calculation yields several equivalent forms for F(z) 



(45) 
(46) 

(47) 
(48) 



F(z) 



E 



(n - 4)fe„ 



£j(„-2)(n-l) „ =5 



7y-E( n " 4 ) fc « 



n — 2 n — 1 



(49) 



OO 



(n - 4)fc n 



x 



_ 5 (n-2)(n-l) 
1 + z + z 2 + • • • + z n - 3 - (n - 2)z"- 2 ] 

(to - 2)fc m+2 - (to - 3)fc m+ i 



(50) 



^ m— 3 



n 



[1 + z + z 2 + + 



(51) 



-W„-2)(n-l) 



n— 5 



[l + 2z + 3z 2 + • • • + z™~ 4 + (n - 2)z 



n — 31 



(52) 



The multiple expressions given above are useful in different situations. The 
first expression is most useful near spatial infinity (where it implies F(z) = 
F(0) + O(z 3 ) and makes manifest the rapid falloff of the remainder of the stress- 
energy once the Hawking flux has been subtracted off). The second, third, and 
fourth expressions are most useful near the horizon, where the explicit (1 — z) 
factor helps to keep the behaviour of this particular term mild. Finally, the 
third form given above can also be cast into the form 



F(z)= Poc (l-z)J2< 



j—max(£-\-l,3) 



(j - 2)k j+2 - (j - 3)k j+ i 



If we work near spatial infinity this is enough to show that 



p(z) 



+ 



p(z) = 0(z 



m = + 



20 P ° 

JL 

20 Pc 
4 ). 



F(0) 
2 

F(0) 



1 - z 



1 - z 



+ 0(z 4 ). 



+ 0{z A ). 



20 Po ° 2 



1 - z 



On the other hand, near the horizon we have 



(53) 



(54) 

(55) 
(56) 
(57) 



p{z) 

r{z) 
p(z) 

m 



= + 



20 1 



2 

F(0) 



1 - z 



0(1). 



20 P °° 2 



1 - z 



+ 0(1). 



= + 



0(1). 



F(0)- 



20 



1 -z 



(58) 

(59) 
(60) 
(61) 
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Note the crucial change in the signs of the energy density and radial tension 
needed to keep the stress-energy regular on the future horizon. 

The second form of F(z) given above can now be substituted into the stress- 
energy to show 



p{z) = 



20 



(n - 4)fc„ 



2(n- 2)(n- 1) 



n— 5 



1 - z 



+ J^Poo Z 2 (l + Z + Z 2 



2 V ( n ~ 4 ) fc " 



z 3 -9z 4 ) + p 00 ^k n z n 



n— 4 



n— 5 



(n- 2)(n- 1) 



[1 



2 + 2 



+ z 



n—3 



(n-2)z r 



(62) 



r(z) = + 



1 \^ 

20 P °° 



(n - 4)fc„ 



-2(n-2)(n-l) 



n— 5 



1 - z 



Poo Z 2 (1 + Z + Z 2 + Z 3 + Z 4 ) + P ooJ2 k nZ U 



n=4 



^ (n- 2)(n- 1) 



,n— 3 



(n-2)z r 



/(*) 



2^) Po ° 



(n - 4)fc„ 



n -j2(n-2)(n-l) 



1 - z 



(63) 
(64) 



From the preceding analysis it is clear that this model satisfies all the known 
properties of the stress-energy tensor in the Unruh vacuum (anomalous trace, 
covariant conservation, asymptotic behaviour both at spatial infinity and the 
horizon), and is the most general form of the stress-energy to do so. Conse- 
quently these equations provide a general formalism for the stress-energy tensor 
in the Unruh vacuum. 



2.7 Model building: the final model 

With these analytic preliminaries out of the way, the final construction of the 
semi-analytic model is quite anti-climactic: Merely truncate the power series for 
p(z) at some finite integer n to obtain a rational polynomial approximation to 
the stress-energy tensor. 

We know that the power series starts off at order z 4 . Further we know that 
the anomalous trace introduces order z 6 terms. So good first guess is to simply 
try the three-term polynomial 
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p(z) = poo z 4 (fc 4 + k 5 z + k 6 z 2 ). 



(65) 



Remarkably, this ansatz is sufficient to fit all the known numeric data within 
expected tolerances. Substituting this ansatz into the formulae above completes 
the model: 



p(z) 



r(z) 



h k 6 \ 

Po 



+ Poo (k 4 z 4 + k 5 z 5 + k 6 z 6 ) 
12 



-Poc z 



(1 + z + z 2 - 3z 3 ) + ^(1 + z + z 2 + z 3 - 4z 4 ) 



+ Poo {k A z i + k 5 z 5 + k 6 z 6 ) 

.,2 



+PooZ 



^-(1 + z + z 2 - 3z 3 ) + ^(1 + z + z 2 + z 3 - 4z 4 ) 
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20 24 20 
A little brute force yields 



1 - z 



(66) 



£ &5 ke \ Z £ 2 , i 2 i 3 i 4 \ 

20-24-20) p ~ — 'To p - z ( 1 + z + z +z +z ) 



• (67) 
(68) 



1 - z 1 



20 24 20 



+ k i z 2 +[-k i + ^-) z 3 



5/C5 3/cfi 



6 \ ^4 



1 - z 



4 2 
20 ~ 24 ~ 20 



9£ _ JM 5 
10 5 ' 



2k, 



+ k A z z + \-k 4 + — ) z 



5 \ ~3 



Performing a least squares fit on the transverse pressure data yields 
p(z) = Poo z 4 (26.5652 - 59.0214z + 38.2068z 2 ). 



(69) 



(70) 



(71) 



The energy density and radial tension are then derived quantities with the values 



p = -— (5.34889 + 26.5652z 2 - 105.260z 3 + 35.0869z 4 + 32.91052 5 ) . 

1 - z 
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(72) 

z 2 

t = (-5.34889 + 26.5652z 2 - 65.9128z 3 + 63.3694z 4 - 13.3241z 5 ). 

1 — z ' 

(73) 

(I am quoting six figure accuracy only for clarity of exposition — we should not 
trust these formulae beyond the 1% level.) 

This data-fitting was performed in the following manner: Through the kind 
offices of Professor Ottewill and Professor Jensen I acquired copies of the nu- 
merical data they used in their 1991 paper jl3| [l4[ |l^|. For spin zero, they 
provided me with the expectation values of the stress-energy tensor in the Un- 
ruh vacuum, which they had calculated by numerically evaluating the difference, 
(U\p(z)\U) — (H \p(z)\H), between the transverse pressures in the Unruh and 
Hartle-Hawking states, and adding this difference to the previously published 
data of Howard |^3| . After checking their data for internal consistency |24| , I 
subtracted Howard's values for (H\p(z)\H) in order to reconstruct these differ- 
ences. I then added these differences to the improved Anderson-Hiscock-Samuel 
values for (H\p(z)\H) |25|, ^6). I discarded the point on the horizon, since the 
Hartle-Hawking data at the horizon is itself an extrapolation. Also, because 
the various data sets are calculated for somewhat different values of r, I had to 
discard several other data points when merging the sets. Finally I added the 
point at spatial infinity because we have exact information at that point. The 
resulting data set contained 26 points, and is summarized in the table given 
herein. (To avoid unnecessary proliferation of numerical data, I am including 
only the bare minimum required to construct the semi-analytic model.) Math- 
ematica was then used to perform a simple least-squares fit, thereby producing 
the coefficients given above. 

By comparing the fitted curve to the original data set, I checked that the 
maximum deviation was 0.8%. Given the expected 1% accuracy of the numeric 
difference data further refinement of this model does not seem currently justi- 
fiable. When graphically plotted the actual numerical data cannot be visually 
distinguished from the fit. 

Once a fit to p{z) has been agreed upon, there is no further flexibility in 
the model, and p(z), t(z), and the flux f(z) [equivalently, the total luminosity 
L] are completely specified. In particular, it is meaningless to try independent 
curve fitting to the energy density or radial tension. 

Indeed, if we obtain a good fit to the transverse pressure, which then does 
not result in a good fit to the energy density or radial tension, this does not mean 
that the fit is bad. On the contrary, since the transverse pressure is used (in both 
the semi-analytic model and the numerical estimates) to calculate the energy 
density and radial tension any discrepancy in these quantities must be ascribed 
to actual error (either numerical roundoff of something more systematic) rather 
than failure of the fit p7] . 

There is an important consistency check that the semi-analytic model should 
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Table II. Numerical data used to create the semi-analytic model. 



r/2M 


Ph 


Pu 


-Ph 


Pu 


1.05 


8.54318 


-4 


426 


4 


11718 


1.10 


7.37818 


-4 


312 


3 


06618 


1.15 


6.56051 


-4 


194 


2 


36651 


1.20 


5.96614 


-4 


076 


1 


89014 


1.25 


5.51801 


-3 


962 


1 


55601 


1.30 


5.16771 


-3 


848 


1 


31971 


1.35 


4.88431 


-3 


740 


1 


14431 


1 /in 
1.4U 


4.04 10\J 


Q 

-o 


Oo4 


1 


01380 


1.45 


4.44508 


-3 


536 





90908 


1.50 


4.26740 


-3 


440 





82740 


1.55 


4.10898 


-3 


350 





75898 


1.60 


3.96552 


-3 


264 





70152 


1.65 


3.83431 


-3 


182 





65231 


1.70 


3.71331 


-3 


108 





60531 


1.75 


3.60103 


-3 


034 





56703 


1.80 


3.49635 


-2 


966 





53035 


1.85 


3.39840 


-2 


900 





49840 


1.90 


3.30648 


-2 


838 





46848 


1.95 


3.22005 


-2 


780 





44005 


2.00 


3.13862 


-2 


724 





41462 


2.10 


2.98924 


-2 


620 





36924 


2.20 


2.85568 


-2 


528 





32768 


2.30 


2.73583 


-2 


444 





29183 


2.40 


2.62795 


-2 


368 





25995 


2.50 


2.53057 


-2 


298 





23257 


Infinity 


1.00000 


-1 


000 





00000 



The Hartle-Hawking data is from Anderson-Hiscock-Samuel Pq, El 
The difference data is inferred from Jensen-McLaughlin-Ottewill |13|, jl 
The final column is the input to the least-squares fit. 



16 




z=2M/r 



Figure 1: This graph shows the transverse pressure in the Unruh vacuum, 
both fitted and numerically estimated, as a function of z = 2M/r over the 
entire available data range. 




z=2M/r 



-10' 



Figure 2: This graph shows the radial tension in the Unruh vacuum, both 
fitted and numerically estimated, as a function of z = 2M/r, over the entire 
available data range. 



17 




Figure 3: This graph shows the radial tension in the Unruh vacuum, both 
fitted and numerically estimated, as a function of z = 2M/r, focussing on small 
z. 




Figure 4: This graph shows the energy density in the Unruh vacuum, both 
fitted and numerically estimated, as a function of z = 2M/r, over the entire 
available data range. 
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density 




z=2M/r 



Figure 5: This graph shows the energy density in the Unruh vacuum, both fitted 
and numerically estimated, as a function of z = 2M/r, focussing on small z. 
Note the discrepancy in the region below z = 0.5, r > AM. (See footnote p7J.) 



pass: Elster has calculated |l6| the total luminosity of a black hole against 
scalar emission using standard techniques |28 which are totally orthogonal to 
the present analysis. At each frequency he calculates a transmission coefficient 
which describes how much of the radiation leaving the event horizon actually 
makes it out to spatial infinity. [In (3+1) dimensions some radiation is always 
back-scattered from the gravitational field.] After multiplying this transmission 
coefficient by a Planckian spectrum and integrating over all frequencies Elster 
reports the total scalar luminosity as 



L Elster = 7.44 x 1(T 5 M 
In the notation of this paper, this translates to 

81 7.44 



El.il < 



Po 



16 6.9941 

On the other hand, the semi-analytic model gives 



5.385 po 



Jo 



96 
20 



59.0214 
24 



38.2068\ 
-20- J P ° 



5.349 Pc 



(74) 



(75) 



(76) 



That two such rather different techniques give total luminosities in such good 
agreement (better that 0.7%) is not only very encouraging, but is in fact a 
powerful consistency check on the numeric data. This completes construction 
of the semi-analytic model. 
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2.8 Matyjasek's analysis 

To compare my model to the model developed by Matyjasek JlTj we should first 
note that Matyjasek does not base his model on the general Christensen- Fulling 
analysis as outlined and developed above, but instead uses a less general ansatz 
based on the Brown-Ottewill approximation |12| , What he calls his N = 6 
ansatz is (as presented) internally inconsistent. It is incapable of simultaneously 
fitting the luminosity data and giving the right asymptotic behaviour at infinity. 

More precisely, if we fix the correct asymptotic behaviour at spatial infinity, 
then the luminosity is not a free parameter in Matyjasek's N = 6 ansatz. His 
ansatz is equivalent to keeping k± a free variable and setting k§ = +6 and ke = 
—9, in which case F (0) = —2/5 and /o = hpoo- (This is close to but not equal to 
the geometric optics approximation f geometric- optics = 81/16 = 5.0625 ^ 5; it 
is certainly not equal to Elster's calculated value fEister — 5.385.) This analysis 
also forces Matyjasek's 62 coefficient to be exactly 62 = 5. (Matyjasek's T is my 
Poo-) With only one free parameter (64) the N — 6 ansatz can at best give only 
qualitative agreement with the numeric data. 

Matyjasek's N = 7 ansatz has enough free parameters to fit both the asymp- 
totic behaviour of the stress-energy and leave the luminosity as a free vari- 
able. This model is equivalent to setting k$ = —9, fcs = 126 — 24/o/poo, and 
fc 4 = b| 2 =i + 24/ ]/p oo -117. So in this model F(0) = -2/ + (96/10)^. Maty- 
jasek's N = 7 model has only two free parameters, /o and p\ z= i, the flux and 
the transverse pressure on the horizon, which he fits to these two single pieces 
of data as calculated by Elster and Jensen-McLaughlin-Ottewill. Matyjasek's 
analysis does not use any of the other numeric data in any quantitative manner. 
With only two free parameters the N — 7 ansatz still gives only qualitative 
agreement with the numeric data. 

In contrast, my model has three completely free parameters, fej, k$, and 
fee- I perform a global unconstrained fit to the totality of the available data, 
and provide a quantitative statement on the quality of the fit (1% accuracy). 
Furthermore, I use Elster's luminosity calculation as a consistency check rather 
than as input. 

3 Energy conditions 

Outside the event horizon, the null energy condition (NEC; (Tf i pki > -k 1 ' > 0?) can 
be completely analyzed by looking at three special cases: outgoing null vectors, 
ingoing null vectors, and transverse null vectors. The NEC reduces to the three 
constraints 

p{r) - r(r) ± 2f > 0? p{r) + p{r) > 0? (77) 
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3.1 Numerical analysis using the semi- analytic model 

Numerically, these conditions are best investigated visually, by inspection of the 
relevant graphs. 

NEC (ingoing) 




Figure 6: p — t + 2/: This graph shows the NEC (evaluated using the semi- 
analytic model) on ingoing radial null vectors as a function of z = 2M/r. Note 
that the curve is everywhere positive, so that the NEC is satisfied on ingoing 
null geodesies. 

The fact that p — r — 2/ is negative over the entire range z <G [0, 1] , is enough 
to imply that at least in the outgoing radial direction the NEC is violated 
throughout the region exterior to the event horizon. Ipso facto all the pointwise 
energy conditions (null, weak, strong, and dominant) are violated throughout 
the entire region outside the event horizon, all the way to spatial infinity. 

The pointwise energy conditions are violated in the sense that at every point 
outside the event horizon there is at least one null or timelike vector violating 
the conditions. It is certainly not true that all null or timelike vectors violate 
the energy conditions, nor is it true that the energy condition violations are 
large. What is true is that the violations are widespread. 

Furthermore, this implies that the averaged null energy condition (ANEC) 
is violated on outgoing radial null geodesies. (More precisely, the one-sided 
ANEC that is cut off at the event horizon is violated). This automatically 
implies violations of (one-sided) averaged weak (AWEC) and averaged strong 
(ASEC) energy conditions. 

One can also immediately see that the null energy condition is violated along 
the unstable photon orbit at r = 3M (z = 2/3). 

We could now add more detailed analyses of exactly which energy condition 
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NEC*(l-z) (outgoing) 



0.2 




Figure 7: p — r — 2/: This graph shows the NEC (evaluated using the semi- 
analytic model and multiplied by (1 — z) to control the singularity at the horizon) 
on outgoing radial null vectors as a function of z = 2M/r. Note that the curve 
is everywhere negative, so that the NEC is violated on ingoing null geodesies. 




Figure 8: p — r — lf: This graph shows the NEC (evaluated using the semi- 
analytic model, multiplied by (1 — z) to control the singularity at the horizon, 
and divided by z 5 to suppress the asymptotic behaviour at infinity) on outgoing 
radial null vectors as a function of z = 2M/r. Note that the curve is everywhere 
negative, so that the NEC is violated on outgoing null geodesies. 
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NEC*(l-z) (transverse) 




Figure 9: p + p: This graph shows the NEC (evaluated using the semi- analytic 
model and multiplied by (1 — z) to control the singularity at the horizon) on 
transverse null vectors as a function of z = 2M/r. Note that the curve changes 
sign near z = 0.6. The NEC is satisfied on transverse null vectors far away from 
the black hole, and violated sufficiently near the black hole. (In particular the 
NEC is violated on the unstable photon orbit at r = 3M.) 

is violated where, along the lines reference ||, but the limited additional insight 
to be gained does not seem to warrant it. 

3.2 Some analytic results 

In the outgoing null direction we have the exact result [cf. eq. (plo|)l 



(U\T**\U) £+£+ = p-T-2f, 



F(z)-F(Q) 



1 4 5 



(78) 



But I have already shown [cf. eq. (^) and the discussion on p. 12] that 
F(z) = F(0) + 0(z 3 ). Therefore, as previously asserted 



(U\T&\U) lpj=0(z 5 



(79) 



In fact, from eq. ( ETj ) we have F(z) = F(0) + z 3 Qi(z), with Qx{z) having a 
power series that starts off with a constant z° term. Thus 



(u\t^\u) em 



~l-z 



Qi(z) 



(80) 
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NEC* (1-z) /z^2 (transverse) 




Figure 10: p+p: This graph shows the NEC (evaluated using the semi- analytic 
model, multiplied by (1 — z) to control the singularity at the horizon, and divided 
by z 2 to suppress the asymptotic behaviour at infinity) on transverse null vectors 
as a function of z — 1M.fr. Note that the curve changes sign near z = 0.6. The 
NEC is satisfied on transverse null vectors far away from the black hole, and 
violated sufficiently near the black hole. (In particular the NEC is violated on 
the unstable photon orbit at r = 3M.) 
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It is this explicit prefactor of z 5 /(l — z) that is responsible for most of the 
structure as seen in figures ^ and 0. 

On the other hand, in the ingoing null direction it is easy to obtain the 
general result 



(U\T^\U)q£- = p-r + 2f, 

i 

5 

£ 

5 



+~ PooZ 2 (l + z + z 2 + z 3 - 4z 4 ) - 2-^—F(z), (81) 

1 — z 



= +~ Poo z 2 (l - z)(l + 2z + 3z 2 + 4z 3 ) - 2^— F(z). 

1 — z 



(82) 



But from eq. fl52] ) we know that F(z) contains an explicit factor of (1 — z) 2 so 
that F(z) = (1 — z) 2 Q2(z) with Q2{z) having a power series that starts off with 
a constant z° term. Thus 



(u\t»"\u) e-^ = z 2 (i- z ) 



Poo(l + 2z + 3z 2 + 4z 3 ) - 2Q 2 (z) 



(83) 



It is this explicit prefactor of z 2 (l — z) that is responsible for most of the structure 
as seen in figure [| 

Finally, in the transverse direction one has the general result 



{U\T^\U)lfl± = p + p, 
= -fo 



±- Pqo z 2 (1 + z + z 2 + z 3 - 9z 4 ) 



1 - z 10 
+2p(z) - 2^—F{z), 



(84) 



T— I -/o + |tPoo (1-z)(1 + 2z + 3z 2 -9z 3 ) 
1 — z 10 



+2(l-js)^-2F(z) 



(85) 



But p(z) = z 4 Q3(z) where Qs{z) has a power series that starts off with a 
constant z° term. Thus 



([/|t^|c/) eUj- = 



/o + ^Poo(l-10z 4 -9z 5 ) 



(1-z) 

+2(l-z)z 2 Q 3 (z)-2F(z) 



(86) 
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It is this explicit prefactor of z 2 /(l — z) that is responsible for most of the 
structure as seen in figures || and 

3.3 Some explicit results for the three-term model 

This generic behaviour can easily be checked analytically for the simple three- 
term model developed in this paper. A little brute force yields 



P - r + 2/ 



PooZ 2 (l - z) 



4 20 24 20 ~~ ' 



2£_ k 5 2k 6 



P-T-2J 



P + P 



z 



3£ k 5 3fc 6 
~5 ~ ~2 ~ T~ 

2fe 



z 2 + 



4^ _4fcg 
5 5 



(87) 



1 - z 

z 1 
1 - z 



20 24 20 



9fc5 5A;6 



z 4 + 



. 5 5 
2k 4 z 2 + (-2k 4 

9^ I2k 6 
10 5~ 



7fc 5 



(89) 



The known analytic behaviour of the prefactors in these expressions serves as a 
consistency check on the numerical analysis used to generate the figures. 



4 Discussion 



In summary, in this paper I have developed a systematic way of building semi- 
analytic models for the stress-energy tensor in Schwarzschild spacetime. For 
the Unruh vacuum I have carried the program forward to the extent of explic- 
itly deriving a three-parameter approximation to the total stress-energy that 
successfully fits all known data to better than 1% accuracy. The model passes 
the consistency test of correctly predicting the luminosity from the fit to the 
transverse pressure. 

In two appendices I sketch how this program can be extended to the Hartle- 
Hawking and Boulware states. 

A central result of this paper is the observation that violations of the energy 
conditions, both pointwise and averaged, are ubiquitous (though small) in the 
Unruh vacuum. This (3 + l)-dimensional result is qualitatively in agreement 
with the (1 + l)-dimensional analytic model considered in |^0|. Furthermore in 
view of the results quoted in ^ ^) we know that this is not a peculiarity of the 
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Unruh vacuum, but that energy condition violations are also widespread in the 
Hartlc-Hawking and Boulware vacuum states. 

Note that I am claiming that the violations are widespread — I am not claim- 
ing they are large. These are intrinsically quantum mechanical effects with the 
typical scale of the effect near the horizon being given by 

- ^ (90) 

Ford and Roman have argued (2^, |3^, |3l| that in many situations the quan- 
tum inequalities may be of more interest than the energy conditions themselves. 
It seems that even if the energy condition violations are widespread, the quan- 
tum inequalities may more stringently constrain the dynamics |3l]]. The semi- 
analytic model developed in this paper may be of some interest in explicitly 
testing the quantum inequalities. (The analysis of this paper is fully consistent 
with qualitative features of the Ford- Roman results of |3(J , but the semi-analytic 
model of this paper will allow one to evaluate the relevant integrals more accu- 
rately. Note that although the coefficients were fitted using data from outside 
the event horizon there is nothing to stop us from taking the resulting model 
and applying it inside the event horizon.) 

Finally, I remind the reader that issues of quantum mechanical violation of 
the energy conditions are of central importance to any attempt at taking the 
classical collapse theorems f22| , the classical topological censorship theorem |B2| , 



or the classical positive mass theorem 33 and trying to see whether or to what 



extent these classical theorems can be generalized into the semiclassical world. 
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A Model building: Hart le— Hawking vacuum 



There is nothing particularly sacred about the Unruh vacuum when it comes 
to model-building of the type discussed in this paper. Similar procedures can 
be applied in the Hartlc-Hawking vacuum as well. This sort of modeling al- 
lows us to get a better handle on the underlying rationale behind the Page 
approximation — this approximation is in some sense (to be described below) 
the minimal form of the stress-energy tensor compatible with regularity on the 
horizon and a "thermal-bath plus curvature-corrections" ansatz at spatial infin- 



A.l General analysis 

For general background information see Jll|. First, since the Hartle-Hawking 
vacuum is to be regular on the both the future horizon and past horizon we 
must have /+ = = /_. (And hence the flux is identically zero: f(z) = 0.) 
At asymptotic spatial infinity we want the stress-energy to look like that of an 
thermal bath of radiation at the Hawking temperature That is: we need to 
have p(z) — > Poo asymptotically as z — > 0, while p(z) — » 3poo, and t(z) — > —poo- 
(For higher-spin fields p(z) — > Sgpoo, t(z) — ► —gpoo, and p{z) — > +<j , f> 00 , with g 
being the appropriate statistical weight.) 

Substituting the result for / back into the general expression for the stress- 
energy, the various components are seen to be 



p{z) = +■$- Poo z 2 (1 + z + z 2 + z 3 - 9z 4 ) + p{z) - - F(z). (91) 

10 1 — z 

£ z 2 
t(z) = -±- Poo z 2 (l + z + z 2 + z 3 + z 4 )+p{z) + - F(z). (92) 



The transverse pressure is still well-behaved all the way down to the horizon, 
so it still makes sense to look for a convergent power-series. But because of the 
known asymptotic behaviour, it is more useful to introduce four new parameters 
(fco,fci,fc2,and &3; not present for the Unruh vacuum) and write: 



ity 



2 



oo 




(93) 



n=0 



A brief calculation yields 




1 - z 



z 




2 



(3&i - 2k 2 ) ln(z) 
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+(1-*) 



V ? ^I"^" [l + z + z 2 + --- + z"- 3 -(n-2)z"- 2 l 
^ (n- 2)(n - 1) L 



.n— 5 



(94) 



Here I have kept the first few terms separate for clarity. Note that &4 does not 
contribute to 

We are really only interested in the quantity 



1 - z 



F{z) = + Poc {-2k (l-z)-3k 1 z + 2k 2 z 2 -k 3 



2 ,^(1-^) 



-(3fci-2fe) 



z 2 ln(z) 



+z 



E l" 
fn - 



1 - z 
(n - 4)fc n 



n— 5 



- (n-2)(n-l) 
+ Po J -2fc (l - z) - 3/ciz + 2/c 2 ^ 2 - fc 3 



[l + z + z 2 + ■ ■ ■ + z n ~ 3 - (n - 2)z"- 2 ] 
z 2 (l-z) 



-(3fci-2fe) 



z 2 ln(z) 
1 - z 



E 

.71=5 



{n - 4)fc K 
(n-2)(n-l) 



+0° 



.n— 5 



(n - 2)(n - 1) 



Which can now be substituted into the stress-energy to show 



(95) 



p(z) 



Poo | + ^z 2 (l + z + z 2 + z 3 -9z 4 ) 
+fc (3 - 2z) + Ak\z - k 2 z 2 + k 3 - 



+(Zh-2k 2 ) Z ^--z> 



1 - z 

(n - 4)fc» 
^(„-2)(n-l) 



(n - 4)fc» 
u ^ ( n-2)(n-l) 



E 



n=4 



E 



[1 + z + z 2 + ■ ■ ■ + z™- 4 - (n - 2)z™- 3 ] L 



(96) 



r(z) 



Pocj - ^ Z 2 (l + Z + Z 2 + Z 3 + Z 4 ) 



29 



-ko(l - 2z) - 2k x z + 3k 2 z 2 - fc 3 



2 , z 2 (l-3z) 



-(3fci -2fc 2 )- 



2 ln(z) 
1 - z 



+ z< 



E 

n— 5 



(n - 4)fc» 
(n - 2)(n - 1) 



n=4 



n— 5 



(n - 4)fc K 
(n-2)(n- 1) 



[l + z + z 2 + ■ ■ ■ + z n - i - (n - 2)z™- 3 ] j. 

(97) 



From the preceding analysis it is clear that this model satisfies all the known 
properties of the stress-energy tensor in the Hartle-Hawking vacuum (anomalous 
trace, covariant conservation, asymptotic behaviour both at spatial infinity and 
the horizon). Consequently these equations provide a general formalism for the 
stress-energy tensor in the Hartle-Hawking vacuum. 

At first glance the presence of the logarithmic terms may be disturbing. Note 
that the logarithms show up only in the combination S(z) = [z 2 ln(z)]/(l — z). 
This combination remains finite both at the horizon (lim 2 _>i S(z) = —1) and 
at spatial infinity (lim^o S(z) — 0), so the logarithmic terms should not be 
excluded a priori. 

There is however a popular ansatz that justifies eliminating the logarithmic 
terms. Near spatial infinity we expect the stress-energy tensor to be that of a 
red-shifted thermal bath of radiation with curvature corrections. Schematically: 

(H\T^\H)(z) = ( H \ T ^\ H )Joo + " curvature -corrections"? (98) 
(1 - z) 2 

Two versions of this ansatz will now be used to more precisely fix the form 
of the stress tensor. 



A. 2 Weak thermal bath ansatz 

Since the curvature is proportional to M/r 3 , we expect the "curvature correc- 
tions" to be (at worst) of order z 3 , in which case we have the conservative 
ansatz: 



(H\T^\H)\ Z = (JTIT^IfOloo (l + 2z + 3z 2 ) + 0(z 3 )? (99) 

If we adopt this ansatz, it enforces very specific choices on the first few coeffi- 
cients. In fact, by looking at p(z) it is trivial to see that fco = lflS and ki = 2g, 
while k 2 = 3g. (Here g is the statistical weight, g = 1 for scalars. Also, re- 
member to use the appropriate value of £ for higher spin.) This is enough to 
make the coefficients in front of the logarithmic terms vanish. It is easy to see 
that the 0(1) pieces of p(z) and t(z) are 3g and —g respectively, while the O(z) 
terms are 6g and — 2g respectively. Additionally, the 0(z 2 ) pieces of p(z) and 
t(z) are proportional to 
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n— 5 



(n-4)fc» 
(n-2)(n- 1) 



3 x 3 5 7 



and 



10 



9g 



2 



(n-4)fc n 
^( n _2)( n -l) 



E 



(100) 



(101) 



respectively. So the conservative ansatz (|99| ) can be simultaneously satisfied for 
all components of the stress-energy provided we pick 



*S = 2 



n— 5 



(n - 4)fc» 
(n - 2)(n - 1) 



(102) 



Note that with this choice of coefficients, not only do we satisfy the ansatz given 
above, but we also guarantee that the stress-energy tensor will agree with Page's 
analytic approximation O, ffl at least to order 0(z 3 ). With this ansatz in place 



P(z) = Po 



,g(l + 2z + 3z 2 ) + k 3 z 3 + k nz n 



n=4 



(103) 



The other components of the stress-energy tensor reduce to 



z 3 £ 

p{z) = pool +3g(l + 2z + 3z 2 ) + k 3 — + -^z 3 (l + z + z 2 -9z 3 ) + Y,k n z n 

{, n=4 

~* t („ ( : 2)(f-" ij (1 + ^ + - + ^ - (n - 2)0 } , 



(104) 

3z 3 £ °° 
t(z) = Poo <| - 3(1 + 2z + 3z 2 ) + fc 3 — - ^ 3 (1 + z + z 2 + z 3 ) + ^ fc„z J 

(n - 4)fc„ 



+* 3 E 1 J is (l + z + z 2 + • • • + z- 4 - (n + 2)z"" 3 ) 
j^(„_2)(n-l) v 7 



(105) 



This completes the weak ansatz. Note that k 3 is not a free variable and that it 
is fc 4 and higher order coefficients that uniquely fix the stress tensor. 
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A. 3 Strong thermal bath ansatz 

A more radical ansatz is to assert that the "curvature-corrections" to the stress- 
energy tensor should be of asymptotic order (curvature) 2 , that is of asymptotic 
order z 6 . In this case we assert 



(H\T^\H)\ Z = (H\T^\H) \ x (1 + 2z + iz 2 + Az 3 + 5z 4 + 6z 5 ) + 0(z 6 )? 

(106) 

We should note in particular that the Page approximation satisfies this more 
radical ansatz. By considering p(z) we automatically deduce k — lg, fci = 2g, 
k 2 = 3g, fc 3 = Ag, fc 4 = 5g, and fc 5 = 6g 7 with k & and higher being left free by 
this ansatz. Instead of using brute force it is useful to first consider a subsidiary 
ansatz for p(z): 

p(z) = Poo 5 (1 + 2z + 3z 2 + Az 3 + 5z 4 + 6z 5 + fc 6 ). (107) 
It is now straightforward to compute 



z 



2 



1 - Z 



F{z) = Poo g (-2 -Az + Az 2 + 2z 3 ) + ^z 2 (l + z + z 2 - 3z 2 ) 

+ ^z 2 (l + z + z 2 + z 3 -4z 3 )]. (108) 



Inserting this into the stress-energy tensor 



p = poo + 2z + 3z 2 + Az 3 + 5z 4 + 6z 5 + fc 6 ) + (+2 + 4z - Az 2 - 2z 3 ) 
- l -z 2 (l + z + z 2 - 3z 2 ) - ^z 2 (l + z + z 2 + z 3 - Az 3 ) 

+ (£M z 2 {1 + z + z 2 + z 3_ Qz 2 ) y (109) 

Collecting terms 

P = + » + «"> 'J- ~ ' ) + ^ ( + 2 + gM^j) 



io 7 V io 



+,« (H^i^iM) j. (110 » 

This can be made to satisfy the strong thermal bath ansatz above if we pick 
- fee - 5 = 100, that is 
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&6 = --105. (Ill) 
9 



Once we do this 



P = Poo.9 \ (1 + 2z + 3z z + Az A + 5z 4 + 6z°) + l± - 105 ) }> . (112) 
p = Poo.g<{3(l + 22 + 3z 2 + 4z 3 + 5z 4 + 6x 5 )+ [ J- - 147 ) }> . (113) 
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r = p 00 .g|-(l + 2z + 3z 2 + 4z 3 + 5z 4 + 6z 5 )+ ^-63^|. (114) 

(The computation for r is completely analogous to that just performed for p.) 
This is exactly Page's analytic approximation [O, H|. 

What does this computation tell us? If we take the strong ansatz, (that 
the stress-energy is a red-shifted thermal bath with curvature corrections of 
order z 6 ), and supplement this with the requirement that there be no terms of 
order higher than z 6 , the we are led uniquely to the Page approximation for the 
stress-energy tensor surrounding a black hole. That is, the Page approximation 
is (in the sense described above) the minimal ansatz compatible with general 
properties of the stress-energy tensor. 

Of course the Page approximation is not exact, and the stress-energy around 
Schwarzschild black holes does have some higher order contributions. Can we 
say anything about 0(z 7 ) terms and higher? Indeed yes, when adding 0(z 7 ) 
terms and higher to p(z) we should be careful to not destroy the thermal bath 
ansatz for lower order terms. Suppose we write 



P{z) =PPage(z) + 5p{z) . (115) 



Then we have 



p(z) = p Pagc {z)-^—5F(z). (116) 

t(z) = TPage {z) + ^—-8F{z). (117) 

So as to not destroy the strong thermal bath ansatz we need both Sp(z) = 0(z 6 ) 
or higher and 8F(z) = 0(z 4 ) or higher. We have already seen that we cannot 
possibly achieve this with any monomial in z. The best we can hope for is to 
find some suitable binomial such as 

Sp n (z)=z n -e n z n+1 . (118) 
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In which case 



5F n (z) = (1 - *){ - ^ 2)(n 4) _ i} [1 + z + z 2 + ■ ■ ■ + z- 3 - (n - 2)z- 2 ] 
- £ " [1 + , + , 2 + • • • + - (n - IK- 1 ] | 

(119) 

So if we pick e„ = [n(n — 4)]/[(n — 2)(n — 3)] there are massive cancellations and 

n — 4 

5F n (z) = z"- 2 (l-z) 2 . (120) 

Ti A 



That is to say 



z 2 „„ , , n - 4 



^„W = sAl-s)- ( 121 ) 



1 — z n — 2 

Consequently the binomials 

W-"- (n-2)(„-3) -" + ' (122) 

(for n > 6) are a useful "basis" for the pieces of the stress-energy tensor that go 
beyond the Page approximation. These basis elements are useful in the sense 
that they do not perturb the lower order pieces of the stress-energy. 

I shall report elsewhere the results of performing such fits to the Anderson- 
Hiscock-Samuel data. 
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B Model building: Boulware vacuum 



Unsurprisingly, a similar analysis can be applied in the Boulware vacuum. 
Again, for general background information see Eq j. First, since there is no 
net flow of radiation in the Boulware vacuum we must have /+ = /o = /_. 
(And hence f(z) — 0). At asymptotic spatial infinity we want the stress-energy 
to be as small as possible. Since if nothing else the anomaly will generate 
terms of order 0(z 6 ), Christensen and Fulling were led to tentatively suggest 
(H\T^\H) = 0(z 6 ) @ page 2096]. 

To at least force the terms of order 0{z 2 ) to vanish we need to enforce 



1 + G(0) + 2/ = 0. (123) 



That is 



A-liP--^. (124) 

(This is formally very similar to what we found for the Unruh vacuum, up to 
a few critical minus signs.) Substituting this result for /o back into the general 
expression for the stress-energy, the various components are seen to be 



p(z) = -^ Poo z e ^—- + 2p(z) + — lG(z)-G(0)}. (125) 

r{z) = +^ Poo ^--^-[G(z)-G(0)}. (126) 

I have actually pulled a minor swindle to get to this point because G(0), 
G(z), and /o are all ill-defined infinite quantities in the Boulware vacuum. This 
happens because p(z) diverges at the horizon, and the integral used to define 
G(z) does not converge. Fortunately this does not matter, since the final ex- 
pression for the stress-energy contains only terms of the type 

G(z) = \G{z) - G(0)] = - £ ( ^ - A) P (z) dz. (127) 

This integral converges provided p(z) is of order 0(z 3 ) or smaller at spatial 
infinity. The divergence at the horizon does not matter because the range of 
integration does not include the horizon. 
We can further define 

F(z) = [F(z) - F(0)\ ^-J*z*(l- z) ^ dz . (128) 

Then, integrating by parts, it is easy to show that 
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G(z) 



1 - z 



p(z)-F(z). 



(129) 



Doing this requires only the mild assumption that p(z) is is of order 0(z 3 ) or 
smaller at spatial infinity, which we already had to assume anyway. 
The components of the stress-energy tensor can now be rewritten as 



/ x C 6 l0-9z 

p(z) = - T ^z^—r 



1 - z 



F(z). 



r(z) 



10 P °° l-z 



-p(z) 



1 



-Hz). 



(130) 
(131) 



The transverse pressure is still well-behaved until you get to the horizon, so it 
still makes sense to look for a power-series; though now we expect this power 
series to have radius of convergence of one, and to diverge at the horizon. Be- 
cause of the minimal 0(z 3 ) asymptotic behaviour I have argued for above, it is 
instructive to write 



p(z) = Po 
A brief calculation yields 



k 3 z 3 + ^2 k nZ r ' 



n=4 



(132) 



F(z) = + Poo {k 3 [z - — 



^2(n-4)k n 



n— 5 



y n—l 



n — 1 n — 2 



Which can now be substituted into the stress-energy to show 



(133) 



P{z) = Po 



r(z) 



£ 6 10-9z 
~ 10 Z l-z 



-h 



2(1 -z) 



-— -^(n-4)fc„ 



n— 5 



z n ~ 1 z n ~ 2 



n — 1 n — 2 



£ z 7 z 3 (4 — 3z) °° 
J> > < +- T7T i — I + fc 3— i — + k nz n 

n=4 



10 l-z 

„2 00 



1 - Z 



+ T — ^(n-4)fc„ 



n— 5 



^n— 1 z"~ 2 



n — 1 n — 2 



(134) 



(135) 



We again see from the preceding analysis that this model satisfies all the known 
properties of the stress-energy tensor in the Boulware vacuum (anomalous trace, 
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covariant conservation, asymptotic behaviour both at spatial infinity and the 
horizon). Consequently these equations provide a general formalism for the 
stress-energy tensor in the Boulware vacuum. 

It is possible to improve the situation by making several simplifying ansatze. 
If we assume that all non-zero components of the stress-energy tensor should 
asymptotically be of the same order, then we must set k 3 = 0. If (following 
Christensen-Fulling) we make the stronger ansatz that that all non-zero com- 
ponents are of order 0(z 6 ) one has k 3 = k 4 = k 5 = 0. In this case it is easy to 
check that 



p(z) = Poc (-£+^k^jz 6 + 0(z 7 ). (136) 

t(z) = Poo 1 k 2 z 6 + 0{z 7 ). (137) 
p{z) = Poc k 2 z 6 + O(z 7 ). (138) 
For k§ = — 16 this is compatible with the Brown-Ottewill analytic approxima- 



tion 1 12 . I have not found any nice principle that would uniquely lead from this 
general decomposition to the full Page-Brown-Ottewill approximation. 

Serious model building is now simply a matter of getting good enough data 
to make reasonable estimates of the various coefficients k n . One might even 
hope that with a little more work it might prove easier to analytically calculate 
these coefficients in preference to the numerically intensive work required to 
calculate the stress-energy components directly. 
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